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I. INTRODUCTION 



Ml properties - nuclear magnetic moments as well as radiative capture cross sections - 
are the fundamental low-energy observables of a few nucleon systems and therefore presents 
ideal laboratory to test effective field theories (EFTs). In this regard, Ml properties have 

nnnnnnQ 

been extensively studied using EFT with huge successes [1|, |2|, y, U, 15|, la, 17|] . One such example 
is the ability to describe cr„p, the capture cross section of the np — *• (i7 process, at threshold 
with 1 % accuracy by applying heavy-baryon chiral perturbation theory (HBChPT) up to 
next-next-next-to the leading order or N^LO [1|. In this work, we extend our up-to N'^LO 
HBChPT description of the Ml properties to A = 3 systems. By taking the magnetic 
moments of and ^He as input to fix the coefficients of the contact-term operators, a 
completely parameter-free theory predictions will be made for the total cross section and 
the photon polarization of the thermal neutron capture process {nd ^H7). We will also 
revisit the theory predictions for the two-body observables: deuteron magnetic moment /i^ 
and (Jnp. 

The purpose of this article is to demonstrate the general tenet of EFTs by studying the Ml 
properties of a few body systems: once the long-range contributions are taken into account 
correctly, EFTs enables accurate and model-independent results, regardless the details of 
the short-range physics. 

We begin with a few comments that are generic to all EFTs. At a certain order in EFTs, 
there appear contact-terms (CTs), which parameterize the high-energy (or short-range) 
physics above the cutoff scale of the theory. The coefficients of CTs - which we refer to 
as low-energy constants (LECs) - are thus sensitive to the short-range physics, and depend 
on the adopted cutoff value and the regularization/renormalization scheme. The values of 
LECs are not fixed by the symmetry alone, and should be determined by either solving the 
underlying theory or by fitting them so as to reproduce selected set of known experimental 
data. Since the former is currently not feasible, the latter remains the only practical option. 
At N^LO, HBChPT Ml currents contain two non-derivative two-nucleon CTs, one in iso- 
vector and the other in iso-scalar channel. These LECs, g^^ and g^si will be determined 
in this work by requiring to reproduce the experimental values of the magnetic moments 
of ^H and '^He. Once these LECs fixed, we are left with no free parameters and can make 
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totally parameter-free theory predictions for the other Ml observables*^. Note that the CT 



contributions have been ignored in Ref. which caused small cutoff-dependence in o"„p. By 
taking into account of the LECs, we will show that cr„p becomes virtually completely cutoff- 
independent. Second comment is about the accuracy of the adopted wave functions at short 
range. The authors of Ref. [81] have developed an approach called EFT* or MEEFT {more 
effective effective field theory) that enables a consistent and systematic EFT calculations on 
top of accurate but phenomenological wave functions. The key observations are following. 
The model-dependence resides mainly in the short-range region of the wave functions. Since 
short-range contributions can be well embodied by local operators at low-energy and since 
EFT has the machinery to contain all the relevant local operators (i.e., CTs) in a consistent 
and systematic manner, the model-dependence due to short-range physics is to be absorbed 
into the renormalization procedure of the LECs. To be more specific, if we adopt other 
wave functions that have different short-range behavior, the values of LECs should also be 
changed so as to reproduce the selected experimental data with the adopted wave functions. 
By performing this procedure, while the values of LECs - which are not physical observables 
- are model-dependent, the resulting net contributions become model- independent. An easy 
and effective way of proving the model-independence in a quantitative fashion might be to 
look at the cutoff-dependence of the results, since the cutoff value is the key parameter that 
characterize the short-range contributions. Such a numerical proof will be taken in this work. 
The third comment is about the long-range contributions. Note that mismatches in the long- 
range contributions cannot be cured by finite set of local operators. The long-range part of 
the transition operator is usually governed by the chiral symmetry, leaving little uncertainty 
there. On the other hand, the long-range part of the wave functions is controlled by the 
effective range parameters (ERPs) such as the nuclear binding energies and the scattering 
lengths. For two-nucleon systems, most of the modern realistic A^A^ potentials reproduce the 
ERPs with a great accuracy. However for nucleon systems with A > 3, situation becomes 
highly non-trivial as many of the available potentials fail to reproduce the relevant ERPs to 
the desired accuracy. 

As we will demonstrate, our results have little cutoff-dependence for all the cases con- 



There are many other ahernatives. For example, one can fix and 54^ from the experimental values of 
(Trip and the deuteron magnetic moment and then make theory predictions on /^('^H) and /i(^He) 
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sidered, which might be interpreted that the short-range physics is well under control. On 
the other hand, the mo del- dependence due to the difference in long-range part of the wave 
functions will cause correlations of the matrix elements with the ERPs. In our work, we 
observe rather a strong model-dependence and demonstrate how it is correlated to the model 
prediction of the triton binding energy B^. It indicates that the mo del- dependence is due 
to the mismatches in the long-range contributions. 

To bypass the difficulty and to get model-independent accurate theory predictions, we 
have explored two different approaches. One is to bring prediction of to its experimental 
value i?3^^ = 8.482 MeV, using the observed correlation curves. The resulting Ml matrix 
elements are found to be model-independent to a good accuracy, and consistent with the 
experimental data. Another way is to adjust the tri-nucleon interactions (TNIs) to meet the 
experimental values of the ERPs. 



II. FORMALISM 



A. Faddeev equations 

During the last few decades several different methods permitting to solve three body 
bound and scattering problem has been developed. In this study we solve Faddeev 9| 
equations (also often called Kowalski-Noyes equations) in configuration space to obtain 3- 
body bound and scattering wave functions. We employ the isospin formalism, i.e., consider 
proton and neutron as two degenerate states of the same particle - nucleon, having the mass 
fixed to h'^/m = 41 AH MeV-fm. Then three Faddeev equations become formally identical, 
having the form 

{E-Hq- V,j) = Vij{<!>jk,i + (1) 

where (ijk) are particle indices, Hq is kinetic energy operator, Vij is two body force between 
particles i and j, ^ij,k is Faddeev component. It is useful to define cyclic (-P^) and anti- 
cyclic (P^) particle permutation operators, which permits to transform Faddeev component 
between two particle bases: P+ = (P^)^^ = P23P12 and P^^ij^k = ^jk,i, while P^^ij^k = 
^ki,j- The wave function in Faddeev formalism is the sum of three Faddeev components, 
which employing permutation operators can be written as: 

^ = (l + P+ + p-)$,^.,. (2) 



Faddeev components, if represented in its proper coordinate basis, have simple structure 
and analytical asymptotic behavior for the short-range potentials. We use relative Jacobi 
coordinates = {vj — ri) and = :^(^fe^~^x^)' whereas Faddeev components we expand 
in bipolar harmonic basis: 

= V ^^^^^hlA (/^ (,^, ) ) . (/ \ ^ I (t^t ) t,) , (3) 

here index a represents all the symmetry allowed combinations of the quantum numbers pre- 
sented in the brackets: and ly are the partial angular momenta associated with respective 
Jacobi coordinates; Si and are the spins and isospins of the individual particles. Func- 
tionals Fa{xk,yk) are called partial Faddeev amplitudes. Three- nucleon system conserves 
its total angular momentum J as well as its projection M, however due to the presence of 
charge dependent terms in nuclear interaction, total isospin of the system T is not conserved. 

Equation ([T]) is not complete, it should be complemented with the appropriate bound- 
ary conditions. Boundary conditions can be written in the Dirichlet form. First Faddeev 
amplitudes, for bound as well as for scattering states, satisfy the regularity conditions: 

F^{0,yk) = F^{xk,0)=0. (4) 

For the bound state problem wave function is compact, therefore the regularity conditions 
can be completed by forcing the amplitudes Fa to vanish at the borders of a hypercube 
[0,X^J X [0,F^ax]: 

-Fa(-^max, Vk) = Fa{Xk, Fmax) =0. (5) 

Finally, we normalize three- nucleon wave function to unity | = 1. 

Faddeev components describing neutron-deuteron scattering, for the energies below the 
break-up threshold, vanish for ^ oo. As ^ oo interaction between particle k and 
cluster ij is negligible and Faddeev components ^jk,i and ^ki,j vanish. Then the component 
^ij^k describes the plane wave of the particle k with respect to the bound particle pair ij: 

lim %fe(xfc,yfc) = 4^y] |{^d(xfc)} . ® {y5„(yfc) ® Sfc} . ) ® {Utj) tk) 

X 2 ^^'""'^^ " ^i^Mi'r^ K ^P'^nd)] , (6) 

where deuteron, being formed from nucleons i and j, has quantum numbers = 1, = 1 
and = and its wave function ipdi'^k) is normalized to unity; p designates the relative 



momentum of incoming neutron, r^d = (\/3/2)7/fc is relative distance between neutron and 
deuteron target, whereas hf^ are the spherical Hankel functions. Expression ([6]) is normalized 
so that nd scattering wave function has unity flux. 

For zero or very low momentum neutrons, as is the case for the thermal neutron capture, 
only relative S-wave amplitudes survives in the asymptote, whereas expression ([6]) simplifies 
to: 

y^^^^'^^jA^k,yk) = -j^^liM^k)}^^® {Yi,Xyk) ^ Sk}j^^ ^^^^ {titj)^Jk^ ^ 



(7) 

Tnd 

where "^^^^and is neutron-deuteron scattering length. For the cases where Urbana type 
three-nucleon interaction (TNI) are included, noting that the TNI among particles ijk can 
be written as sum of three terms Vijk = Vij + VJk + ^ki^ modify the Faddeev equation ([1]) 
into: 



B. Electromagnetic current 

For three-body system one has three 1-body currents associated with each particle and 
three 2-body currents associated with each pair of particles. Thus 

3 

Jem = {Jib + "^25 )• (9) 

i=i,i^{j<k) 

Since the wave-functions |\E') in isospin formalism is fully antisymmetric, the matrix element 
of the current operators can be written as 

3 

i=l,i^(j<k) 

We use the electromagnetic current operators derived from HBchPT, which contain the 
nucleons and pions as pertinent degrees of freedom with all other massive fields integrated 
out. In HBchPT the electromagnetic currents and Ml operator are expanded systematically 
with increasing powers of Q/A^ , where Q stands for the typical momentum scale of the 
process and/or the pion mass, and ~ 47r/^ ~ m ~ 1 GeV is the chiral scale, /tt ~ 92.4 
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(a) (b) (c) (d) (e) 

FIG. 1: Tree diagrams for the electromagnetic current operators. Soft one-pion-exchange, the sum 
of the "seagull" (a) and the "pion-pole" (6) diagrams contribute to the Jij^. Diagrams (c) — (e) 
contribute to the Ji^c at N'^LO. The dot represents the vertex corrections coming from NLO or 
N^LO lagrangian. 



MeV is the pion decay constant, and m is the nucleon mass. We remark that, while the 
nucleon momentum p is of order of Q, its energy (~ ^) is of order of Q'^/m, and consequently 
the four-momentum of the emitted photon g'^ = {u, q) with \q\ = uo also is counted as 
0{Cp' jm). Current operators are obtained up to N'^LO. Note that three-body currents are 
N^LO or higher order, and do not enter in our work. 

Let us list the relevant current operators. The explicit form of magnetic moment operators 
can be found in the ref. \{\. The one-body current including the relativistic corrections reads 



m \ 2m 
uj{2jji - Q 



H iq X cr, 

2m 



w{2fj,i - Qi 
8m^ 



{2ipi X cTi) - 

l^'i Qi 



f^i Qi 



Qi -2 



2m^ 
{Aiq X piai ■ pi 



16m^ 



2qq ■ pi) + (higher orders) 



where Qi and /ij represent the charge and magnetic moment of i-th nucleon, and p = V 
—i V) should be understood to act only on the nuclear wave functions. 

Corrections to the IB operator are due to the meson-exchange currents (MECs). Up to 
N^LO, as mentioned, only two-body (2B) contributions enter. It is to be emphasized that 
MECs derived in EFT are meaningful only up to a certain momentum scale characterized 
by the cutoff A. In our work, we adopt a Gaussian regulator in performing the Fourier 



'^'^ It is worth mentioning that there is a different power counting scheme where the nucleon mass is regarded 
as heavier than the chiral scale, m ~ A^/Q, see refs. 10|,lll| for details. However, the use of this alternative 



counting scheme would not affect the results to be reported in this work since the difference between the 
two counting schemes would appear only at higher orders than explicitly considered here (N'^LO). 
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transformation of the MECs from momentum space to coordinate space jsj. It is to be 
noted that the contributions due to high momentum exchanges (above the cutoff scale) are 
not simply ignored but, as we will discuss later, they are accounted for by the renormalization 
of the contact-term coefficients. 

We decompose the two-body current into the soft-one-pion-exchange (Ivr), vertex correc- 
tions to the one-pion exchange (IvrC), the two-pion-exchanges (27r), and the contact-term 
(CT) contributions. 



J, 



ijk) 



jii^) _|_ TU"-; _|_ Tun-; i j 



(jk) 



2B 



jiik) 
CT ■ 



(12) 



It is noteworthy that there can be additional corrections to the 2-body current coming from 
the so-called fixed term. The fixed-term contributions represent vertex corrections to the 
soft-one-pion-exchange and fixed completely by Lorentz covariance. Because the fixed terms 
make the calculation highly involved, but only give very small contributions in Ml operator 
according to our previous study [71] , we neglected the fixed term contributions in the present 
work. 

The soft-one-pion exchange current is NLO and can be written in terms of Rjk = 
\{rj + rk),r = Vj - r^, f = r/|r|, Sjk = Scr^ ■ f cr^ • r - (Tj ■ CTfc, 



9\ 



(tj X fkYr 



2 V'J 



(x) 



Sjk I g^/lAl'^J -1/oaI 



WlylAr) 



where 



T, 



(0) 



5,ifc 

= ^kY 



f f ■ {(Tj CTfc) - -{(Tj CTfc) 



= ±, X , and the regulated delta and Yukawa functions are defined as 



?/oa(0 
yiAir) 



(27r)3 
d^k 

(2^ 
d 



^2 _|_ 

r did 



ml dr r dr 



yoA[r) 



(13) 



(14) 
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(e) 
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(d) 



"1 r 

I I 
I I 



(h) 



FIG. 2: Diagrams which contribute to J2tt (a)-(i) and Jct (j) at N^LO. 



The one-loop vertex correction to the one-pion exchange has been investigated in detail 
inrefs. 

T^8j7'-A[3-' 5 i/oA 2-'t y2Aj 
lij2NwzTl ■ T2 [(CTI + (T2)yof, + (3ff ■ (cTi + CTs) - (cTi + cra))?/^^] |, (15) 



The values of the LECs {c^, c/\, N^z) should in principle be fixed either by solving the un- 
derlying theory, QCD, or by fitting to suitable experimental observables. Since this has not 
yet been done, we adopt here the estimates given in refs. based on the resonance sat- 
uration assumption and the Wess-Zumino action, (c^^, ca, iVvi/z) — (0.1021,0.1667,0.02395). 

The two-pion exchange diagrams give rise to 

= l£^ 0^ f^^^^^ " ^'^^'^'^ + ^^^^ " )^T(r)] - (r, X r,Yf^^L,{r)^ 



(16) 



where 



Ls{r) = -^X^K, + ^-f{-2K, + AK^+r^Ko + 2r^K,) 
6 dr 6 ar ar 

Loir) = 2K2 + g\{8K2 + 2Ki + 2Ko) - g\{16K2 + hK^ + hK^) + gl-^i^K^), (17) 



and the loop functions K^s are defined in refs. 



fla. 



Finally, contact-term contributions have the form 

= e-'^-^-^g X [g,s{<T, + cr,) + g,yT^^^]5j,{T) (18) 

where §4,3 = rripg^ and 04^ = —mp{G^ + \E^'^). We remark that, three contact terms were 
introduced in refs. Q], whose coefficients are denoted as (74, and E^'^. However, due 
to Fermi-Dirac statistics, only two of them are independent, and consistent with Eq. (|T8j) . 
A similar reduction has been noticedjor the Gamow- Teller operator, where only one linear 
combination of two CTs is required 



C. n-d radiative capture 



In center of mass frame, each currents can be written in the form of 



-iq x I 



iq^^ + jcj 



(19) 



where x is Vi for and Rjk for 



iJk) 



To calculate neutron radiative capture observables we will use multipole expansion. First 
we introduce a shorthand notations for the multipoles: 



FjM{r) = ij{qr)YjM{r) 



Fjlir) 



(20) 



where jiil'^) is the spherical Bessel function; Yjm and 3^j£i are spherical and vector-spherical 
harmonics respectively; r is a vector describing the particle (nucleon or meson), which 
interacts with EM field. Then the electric and magnetic multipoles read 



MjM = Ff^{r)-j, + zq 



J+1 
2J+1 



1/2 



J 



2J+1 



1/2 



1/11 



S 



JM 



/ + 1 \ 



J 



1/2 



FfU^ir) 



■jc + qFZir)-j,. (21) 



2J+ 1/ ' V2J + 1, 

With the explicit expressions of the Ff'[{f), Ml multipoles can also be written as 

jo{qr)j,, - ^j2{qr) {jV -r{r- j^)} 




MiM = i\ ^jiiqr) [r X jj + 

oTT Vovr 



In terms of the reduced matrix elements (RMEs 



qfj'N 



A7r(%n\A:jM\\^tt/, 



(22) 



(23) 
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where Xjm = {-Mjm, £jm), the total nd capture cross section is given by 



2 a f he \ f q \ ^ 



{Vrel/C) \2mcy \hc) U ^ 



2 



+ 



9 

1+ ^^A tH _ 3 



(24) 



Thermal neutron capture proceeds only from doublet Jj^ = | and quartet Jj^ = | 
scattering states, since only these two states comprise nd S'-wave asymptote and thus domi- 
nate low energy scattering. Since final state (the triton) is J^^ = therefore only magnetic 

- — 1 1 - — 1 1 

dipole transition elements m2 = Aif"^ , = Ai'("^ and electric quadrupole transition ele- 

~3 1 

ment 64 = 82'^ do not vanish. Notice that magnetic dipole moments are purely imaginary, 
while electric quadrupole moment is real. 

Experimentally, in addition to capture cross section, photon polarization parameter Rc 
can also be measured. This parameter is given by [3] 

1 



Rr 

3 



7 

2 


m4 


^ + VSRe [mam*] + 


5 
2 


64 


1^ + y24Im [mse^] - VSIm [m4e^] 




ni2 


' + I 


m4 




64 


2 



1 



(25) 



Calculations using expression (fTOl) are numerically stable for all the two- and one-body 
current terms except the ones entering into impulse approximation of Ml operator. This 
issue has been observed and the special numerical procedure developed in reference |l^ , we 
have successfully followed it. 

III. RESULTS 

A. Binding energies and scattering lengths 

In this work we have performed rigorous calculations for several qualitatively different 
realistic nuclear Hamiltonians, which are based on NN potentials defined both in config- 
uration and momentum spaces. Argonne Avl8 |l^ is an accurate local A^A^ potential in 
configuration space. Semi-realistic configuration space potential INOY has been recently de- 
rived by Doleschall [IQ], which can describe binding energies of three-nucleon systems with 



only two-nucleon forces. ISUJ [17|- a recent revision of INOY - further improves description 
of np and pp data and at low energies provides solution for the long standing "Ay puzzle" 
of N-d scattering. We have also tested some chiral N'^LO potentials defined in momentum 



space: Idaho group potential 18| (referred to as I-N3L0), and three different parameteri- 
zations of chiral N^LO potential of Bonn-Bochum group y^]. In particular Bonn-Bochum 
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TABLE I: Values for np singlet scattering length and deuteron binding energy obtained using 
Bonn-Bochum group potentials. 



Model 


^a„p (fm) 


Bh, (MeV) 


B1-N3L0 


-23.60 


2.215 


B2-N3LO 


-23.72 


2.218 


B3-N3LO 


-23.64 


2.220 


Exp.: 


-23.74 


2.225 



group potentials parameterized with set of cut-off values |A,A| = {450, 500} , {450, 700} 
and {600, 700} MeV have been used and are referred to as B1-N3L0, B2-N3LO and B3- 
N3L0 respectively. 

All the NN potentials mentioned above describe the A^A^ data quite accurately. And all 
but Bonn-Bochum group potentials reproduce experimental deuteron binding energy i?^ and 
the singlet np scattering length ^a„p with at least four significative digit accuracy. Values of 
these observables obtained using Bonn-Bochum group potentials are summarized in Table [B 

Our three-body calculations have been carried out considering isospin breaking effects, 

which allow admixture of total isospin T = 3/2 in the wave functions. The Argonne UIX 

□ 

three-nucleon interaction [19| also has been taken into account in the combination with Avl8 
A^A^ potential. 

The relevant properties of three-body systems obtained with the adopted models are 
summarized in Table [TH These values are in perfect agreement with ones obtained by the 



other groups [l^,(2l|-(24]- In 25| we have already published three-nucleon properties for 
INOY and Avl8 models, the small difference in fourth digit of those results compared with 
current ones is due to the small admixture of isospin T = 3/2 states. One should note that 
only INOY, ISUJ and Avl8+UIX models reproduce experimental three-nucleon binding 
energies as well as neutron-deuteron doublet (J = |) scattering length accurately. Chiral 
potentials at N'^LO comprise already two irreducible three-nucleon interaction diagrams 
with contact terms. The strength of these contact terms may be a priori adjusted so as 
to reproduce three-nucleon binding energy and scattering length j26|]. In this work however 
only two-nucleon interaction part of N^LO models was considered. 
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TABLE II: Three-nucleon properties as calculated with different realistic Hamiltonians. They 
contain: nd doublet i^&nd) and quartet [^and) scattering lengths in fm; bound state properties 
comprising binding energy (BE), average kinetic energy ((T)) in MeV's and rms radius rj-ms = \/v^ 
in fm. These values are compared to other theoretical calculations and experimental results. 

nd 3jj Sjje 



Hamiltonian 


Ref. 






BE 


(T) 


^rms 


BE 


(T) 


^rms 


Avl8 


this work 


1.266 


6.331 


7.623 


46.71 


1.769 


6.925 


45.67 


1.810 




[21, 22] 


1.248 


6.346 


7.623(2) 






6.924(1) 






AV18+UIX 


this work 


0.598 


6.331 


8.483 


51.29 


1.683 


7.753 


50.23 


1.716 




[21, 22] 


0.578 


6.347 


8.478(2) 






7.748(2) 






INOY 


this work 


0.551 


6.331 


8.483 


33.00 


1.666 


7.720 


32.22 


1.704 




[16] 






8.482 






7.718 






ISUJ 


this work 


0.523 


6.330 


8.484 


32.95 


1.667 










m 






8.482 






7.718 






i-JN3LO 


this work 


1.101 


6.337 


7.852 


34.54 


1.760 


7.159 


33.83 


1.797 




[23] 






7.854 












B1-N3L0 


this work 


1.263 


6.334 


7.636 


33.60 


1.816 


6.904 


32.79 


1.860 




m 






7.64 












B2-N3LO 


this work 


1.024 


6.339 


7.930 


31.70 


1.777 


7.210 


31.01 


1.815 




m 






7.97 












B3-N3LO 


this work 


1.781 


6.329 


7.079 


47.25 


1.863 


6.403 


46.17 


1.909 




m 






7.09 












Exp: 




0.65±0.04 [20] 6.35±0.02 [20] 


8.482 






7.718 







B. Magnetic moments and thermal neutron capture 



In Table ITIl we present Ml RMEs obtained for INOY Hamiltonian with A=700 MeV, 
listing the contributions from each chiral order. Note that the one-body contribution of the 
iso-scalar Ml RME, m2, is strongly suppressed due to the pseudo-orthogonality between 
initial and final wave functions. The chiral convergence is however not much illuminating, 
i.e., N^LO contributions appear about the same size of NLO. This behavior is mainly due 
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TABLE III: Matrix elements calculated for magnetic moments and thermal neutron capture. These 
results are obtained using INOY Hamiltonian with A=700 MeV. 







m('h) 


/i(3He) 








C 2 ' 2 
^2 


LO: IB 


0.8593 


2.6567 


-1.8100 


395.5000 


-13.6196 


13.1149 


-0.0741 


N^LO: IB 


-0.0057 


-0.0199 


0.0080 


-0.1653 


0.4106 


0.1048 


0.0032 


NLO: Ivr 


0.0000 


0.1515 


-0.1501 


7.0970 


-2.5712 


-0.4289 


0.1562 


N^LO: IvrC 


-0.0029 


0.0839 


-0.0926 


3.1860 


-2.7674 


-0.3465 


0.0000 


N^LO: 27r 


0.0000 


0.0374 


-0.0362 


1.1290 


-1.2504 


-0.1223 


-0.0019 


g45 


0.0338 


0.0457 


0.0449 


0.0000 


-0.9855 


0.2647 


0.0000 


g4V 


0.0000 


0.0733 


-0.0712 


2.3130 


-2.5179 


-0.2267 


0.0000 



TABLE IV: Values of contact term coefficients g4s and g4y , which are obtained by fitting magnetic 
moments of triton and ^He, for INOY Hamiltonian. 

A (MeV) g4s g4y 

500 0.2747 1.8746 

700 0.2313 0.8021 

900 0.1997 0.4613 



to the accidental cancelation between two NLO contributions, the seagull and pion-pole 
diagrams [s^. 

As explained, Ml currents contain two non-derivative contact-terms at N'^LO. Since the 
coefficients of them, g^s and g^v, cannot be determined from the underlaying theory yet, 
we fit these constants by requiring that magnetic moments of and ^He are correctly 
reproduced. The resulting values obtained with INOY potential are given in Table IIV[ 
We remark that (745 and g^v depend on cutoff A as well as on particular choice of nuclear 
Hamiltonian. 

Table |V] shows the cutoff dependence of our results. One-body contributions are cutoff 
independent by their construction. NLO results bring sizable cutoff-dependence, indicating 
that some important pieces are omitted at this level. As is indicated in the table, going 
N^LO but without taking the CTs does not help in resolving the situation. It is only after 
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TABLE V: Dependence of Ml observables for two and three-nucleon systems on cutoff value A. 
These results are obtained using INOY Hamiltonian. 











LO 






A (MeV) 


/i(2H) 


M(3H) 


/i(3He) 


anp (mb) 


and (mb) 


Rc 


_ 


0.8593 


2.657 


-1.810 


309.7 


0.2785 


-0.2369 










NLO 






A (MeV) 


/i(2H) 


^(3H) 


/i(3He) 


anp (mb) 


and (mb) 


Rc 


500 


0.8593 


2.760 


-1.913 


318.7 


0.2972 


-0.3026 


700 


0.8593 


2.808 


-1.960 


320.9 


0.3296 


-0.3538 


900 


0.8593 


2.829 


-1.980 


321.9 


0.3480 


-0.3753 








N^LO without contact term 






A (MeV) 


/i(2H) 




/i(3He) 


anp (mb) 


and (mb) 


Rc 


500 


0.8499 


2.836 


-2.011 


324.1 


0.3612 


-0.3896 


700 


0.8507 


2.910 


-2.081 


327.5 


0.4237 


-0.4366 


900 


0.8515 


2.937 


-2.105 


328.8 


0.4526 


-0.4504 










N^LO 






A (MeV) 


/.(2H) 


/z(3H) 


/i(3He) 


anp (mb) 


and (mb) 


Rc 


500 


0.8584 


2.9790 


-2.1276 


330.9 


0.5012 


-0.4659 


700 


0.8585 


2.9790 


-2.1276 


330.5 


0.4946 


-0.4649 


900 


0.8583 


2.9790 


-2.1276 


330.4 


0.4959 


-0.4650 


Exp.: 


0.8574 


2.9790 


-2.1276 


332.6 ±0.7 


0.508±0.015 


-0.420±0.030 



the CTs taken into account that the results become almost independent of the cutoff, which 
implies that the CTs are quite effective in renormalizing away the details residing in the 
short-range region. 

Results with varying model Hamiltonian are given in Table IVIt with some relevant low- 
energy properties of the potentials. From the table, one observes that /i(^H) and anp are 
rather insensitive, Rc is moderately sensitive and the nd capture cross section, and, is highly 
sensitive on the model Hamiltonian. To understand the sensitivity, let us consider the 
model-dependence of the effective-range parameters (ERPs), which govern the long-range 
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part of the RMEs. The most important ERPs are the binding energies and the scattering 
lengths, which should strongly influence Ml RMEs through the coupling of the long-range 
parts of the three-nucleon wave functions. And indeed, as shown in Fig. [3], the Ml RMEs are 
strongly correlated with the triton binding energy B-^. The correlation is found to be almost 
perfect for m4, while with some fluctuation for These behavior can be explained with 
simple arguments: let us first concentrate on the quartet RME, m^. In spin-quartet states, 
Pauli principle inhibits three nucleons from gathering altogether, and thus observables are 
insensitive to short-range part of three-nucleon interaction. As a result, the nd quartet 

TABLE VI: Predictions for the deuterons magnetic moment and the observables of the thermal 
neutron capture on protons and deuterons. These calculations have been realized by fixing contact 
terms of the meson exchange current in order to reproduce magnetic moments of the triton and '^He. 
These values turns to be insensitive to the cut-off parameter in the interval A = (500, 900) MeV; if 
however variation was larger than one affecting the fourth significant digit it is given in parentheses. 
The constructed AV18+UIX* model gives 2a„d=0.623 fm, ^a„d=6.331 fm and BE(3He)=7.718 
MeV; the I-N3L0+UIX** results are 2a„d=0.634 fm, '^a„d=6.339 fm and BE(3He)=7.737 MeV. 
Both these models are adjusted to reproduce experimental triton binding energy of BE('^H)=8.482 
MeV 



Model 




anp (mb) 


and (mb) 


Rc 


AV18 


0.8575 


331.9(1) 


0.680(3) 


-0.435 


AV18+UIX 


0.8604 


330.6(2) 


0.478(3) 


-0.458 


INOY 


0.8585 


330.6(2) 


0.498(3) 


-0.465 


ISUJ 


0.8585 


331.1(2) 


0.501(2) 


-0.466 


I-N3L0 


0.8574 


330.4(3) 


0.626(2) 


-0.441 


B1-N3L0 


0.8577 


328.7(6) 


0.688(4) 


-0.438(1) 


B2-N3LO 


0.8588 


331.0(4) 


0.609(4) 


-0.448(1) 


B3-N3LO 


0.8549 


330.9(7) 


0.879(8) 


-0.411(2) 


AV18+UIX* 


0.8614(1) 


330.9(3) 


0.476(2) 


-0.457(1) 


I-N3L0+UIX** 


0.8590(1) 


329.7(3) 


0.477(3) 


-0.468(1) 


Exp.: 


0.8574 


332.6 ±0.7 [27] 


0.508±0.015 [28] 


-0.420±0.030 [29] 
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scattering length '^and has httle model-dependence; all the models considered here reproduce 
^a„d in excellent agreement with the experimental data. This explains the perfect correlation 
of 7714 with i?3. On the contrary, spin-doublet states are free from the exclusion principle 
and sensitive to the short-range three-nucleon interaction. This makes the scattering length 
'^ttnd largely model-dependent, see Table IVTl and we might expect that m2 depends not only 
on B3 but also on "^and- However ^a^^ and -B3 are correlated, which is known in terms of the 
Phillips line (3l|. The correlation of "^and with B3 is not perfect, showing small deviations 
from the Phillips line. These arguments are in good accordance with what we observe in 
Fig. [3], which shows the correlation of m2 with respect to with some scatters. 

For noble two-body processes, effective range expansion technique often allows us even 



32| for the Gamow- 



algebraic relation of the RMEs in terms of ERPs, see, for example, Refs. 
Teller matrix element of the p + p-^d + e'^ + u process. The problem at hand is however too 
complicate to allow such a mathematical rigor, and we will limit ourselves to an empirical 
curve fitting. We take the trial function as 

m» ^ (26) 

with 

MBs) = m° + t„ m/Brn" - 1] , (27) 

where the superscript i is the model index; that is, mi*^ {71 = 2,4) and -83*'* stands for the 
RMEs and BE obtained with the i-th model potential, respectively. Varying the value 
of z/, values of and t„ are searched by a chi-square fit. The resulting chi-square is found 
to be parabola shape with minimum at around u = —2.5. The solution with u = —2.5 is 

02(^3) = (-21.87 ± 0.24) - 10.76 [{B-,/ B^'"'')-^-^ - l] , 

04(^3) = (12.24 ± 0.05) + 11.35 [{B^/ B^""")-^-^ - l] . (28) 

The solution is drawn in solid line in the figure. The above curve fitting procedure turns 
out to be quite robust. For example, the curves and the values of m° with u = —1.5 
are almost the same as those with u = —2.5. Even if we try a simple-minded linear fit, 
z/ = 1, we have 02(^3) = -21.73 + 33.65(^3 - 1) and ^4(^3) = 12.14 - 35.69(a;3 - 1), 
X3 = Eyls/E^^. Thus the values of (t>n{B^^^) = m° are quite insensitive to the fitting 
parameter u. Furthermore, with the resulting values of (j)n{B^^^) = m^, we have Rc=- 
0.462±0.03 and (T„d=0.490±0.008 mb, which are close to the experimental data. Therefore 
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one can conclude that the observed strong model-dependence in Ml properties of three- 
body systems can be traced to the different model predictions of B3, and that, once we have 
correct -B3, the theory predictions should be very close to the experimental data with little 
model-dependence. 

We have also tried to adjust the nuclear potentials to have correct ERPs. As mentioned, 
-B3 and "^and are the relevant ERPs. But since the two ERPs are strongly correlated to 
each other, simultaneous reproduction of both is rather tricky. This correlation in particu- 
lar strong due to on-shell NN interaction part, nevertheless three-nucleon interaction can 
break it. Note that UIX TNI potential consists of two terms. In our calculation, we have 
readjusted the parameters of those terms to reproduce S3 and "^and simultaneously with 
the Avl8 and I-N3L0 A^A^ potential. We refer respectively the resulting Hamiltonian as 
AV18+UIX* and I-N3L0+UIX**. In addition some charge dependence has been added to 
UIX*, permitting Avl8+UIX* to reproduce also ^He binding energy. The corresponding 
results are given in the bottom lines of the Table IVIi The most important observation to 
be made is that, while the results of Avl8, Avl8+UIX and I-N3L0 differ dramatically, the 
modified Hamiltonians Avl8+UIX* and I-N3L0-I-UIX** give us almost identical results, 
which confirms the argument that our theory predictions are model-independent once the 
ERPs are correctly encoded. The resulting and and Rc are close to the experimental data, 
but with discrepancy of about two sigmas of the data. 

Before closing this section, we would like to make comparison with other calculations 



for the processes considered in this paper. Viviani et. al. [33|, |34] has calculated the Ml 
properties of A = 2, 3 systems with the currents deduced from the adopted nuclear potentials 
using gauge invariance, adding model-dependent pieces for those part that are not fixed by 
the gauge symmetry alone. Their results have some variations depending on the adopted 
potentials |33| and the details of the treatment of the currents. Without model-dependent 
current part capture cross section is underestimated and = (0.418 ~ 0.462) mb, nevertheless 
one gets Rc = —(0.429 ~ 0.446) quite close to experimental value [3J]- Mo del- dependent 
currents enable to reproduce experimental cross section, however the photon polarization 
parameter R^ = —(0.469) becomes larger than the experimental data. Currents related with 
three-nucleon force further increase capture cross section and photon polarization parameter. 



A very similar to ours calculation has been recently performed by Pastore et. al. 35], in which 
electromagnetic current operators have been obtained up-to N'^LO within EFT framework. 
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A-isobar as well as pions and nucleons are treated as pertinent degrees of freedom. And 
they have applied the currents up-to N^LO to A = 2 and A = 3 systems. To this order, 
the CT terms - that play a crucial role in removing the model-dependence at short-range 
physics - do not appear, and they have observed a large cutoff dependence with a substantial 
under-predictions for and and Rc, CTnd = (0.450 ~ 0.315) mb and Rc = —(0.437 ~ 0.331) for 
the momentum cutoff A = (500 ~ 800) MeV. We also acknowledge that, using the so-called 
pionless EFT approach, Sadeghi et. al. ^ have performed up-to N^LO (in their counting 
scheme) calculation for the and and Rc, achieving a perfect agreement with the data. In 
their calculations, the np cross section as well as the nd scattering lengths and the binding 
energies (of A = 2 and A = 3 systems) are taken as inputs needed to fix their parameters, 
the magnetic moments have not been considered. Since magnetic moments are sensitive to 
the D-wave components of the wave functions, it may not be trivial to have accurate theory 
predictions for the magnetic moments using the pionless EFT. A further study in this issue 
will be extremely interesting. 

IV. DISCUSSIONS 

The most natural candidate for the remaining small discrepancy might be the three- 
body current contributions, which are N^LO. It is not difficult to notice that the leading 
three-body contributions are suppressed for both Ml currents and the nuclear potentials, 
for exactly the same reason. Furthermore, the soft one-pion-exchange appears as the leading 
two-body contributions for both of them. Thus we expect that the ratio of the three-body 
contribution to the two-body contribution is the same order for the Ml RMEs and the 
nuclear potentials. 



(0.05 ~ 0.1). (29) 



M2B {V)2B 

Since the TNIs play a crucial role in reproducing the ERPs of three-body systems accurately, 
we may naively guess that the same will also be true for the relation between three-body 
currents and the Ml properties. More quantitatively, eg. (1291) with Table III tells us that 
the three-body current contribution will be about (2 ~ 4) % for 1712 and 777-4, which is just 
the needed size to remove the discrepancy of and and Rc. The same has been demonstrated 



by Viviani et al. 



34| . where 3-nucleon currents have let to increase neutron thermal cap- 
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FIG. 3: Radiative capture of thermal neutron by deutron: correlation of Ml doublet and quartet 
RME's with triton binding energy. 

ture cross section by 0.033 mb. In this regard, taking into account the three-body current 
contribution - while ignoring other pieces of N^LO for simplicity - might be extremely 
interesting. 

V. CONCLUSION 

In this paper Ml properties, comprising magnetic moments and radiative capture of ther- 
mal neutron observables, are studied in two- and three-nucleon systems. We utilize meson 
exchange current derived up to N^LO using heavy baryon chiral perturbation theory a la 
Weinberg. At N3L0, two unknown parameters, g^g and g^y, enter as the coefficients of 
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contact terms. Following the MEEFT strategy, we have fixed them by imposing the renor- 
malization condition that the magnetic moments of tritium and ^He are reproduced. Then 
we analyze the predictions for other Ml properties: magnetic moment of deuteron, as well 
as observables of the thermal neutron capture on proton and deuteron. Analysis comprise 
several qualitatively different realistic nuclear Hamiltonians, which allows us to judge on 
the model dependence of our results. We obtain stable, cut-off independent results, which 
reconfirms efficiency of MEEFT procedure. Model predictions for two-body observables 
(deuteron magnetic moment and thermal np capture cross section) scatter closely around 
the experimentally measured values. 

Radiative capture cross section of thermal neutron on deuterons varies a lot from one 
Hamiltonian to the other. We have demonstrated that this variation is mostly due to the 
correlation of the capture cross section with a model predicted three-nucleon binding en- 
ergy. By fixing three-nucleon binding energy to the experimental value one can reduce model 
dependence below 2% level and obtain model-independent predictions for thermal capture 
cross section and = 0.490 ±0.008 mb and photon polarization parameter Rc = —0.462 ±0.03. 
Within these model-dependent error bars capture cross section agrees with experimentally 
measured value 0.508 ± 0.015 mb [28]. However photon polarization parameter is ob- 
tained slightly too large, like in other studies based on realistic nuclear Hamiltonians and 



currents [3j]. The remaining discrepancy is comparable in size with higher order terms of the 
EFT, which have been neglected here. We believe that in particular three-nucleon currents, 
which first appear at N^LO in our power counting scheme, should be important. 
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